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Abstract 

Nonhermitean hamiltonians of convection-diffusion type occur in the descrip- 
tion of vortex motion in the presence of a tilted magnetic field as well as in 
models of driven population dynamics. We study such hamiltonians in the case 
of rectangular barriers of variable size. We determine Lyapunov exponent and 
wavenumber of the eigenfunctions within an adiabatic approach, allowing to re- 
duce the original d~2 phase space to a d=l attractor. 
PACS numbers:05.70.Ln,72.15Rn,74.60.Ge 

1 Introduction 

Quantum mechanics in an imaginary magnetic field (IMF) has been studied by var- 
ious authors, after Hatano and Nelson pointed out its relevance in describing the 
competition between drift and pinning for a vortex in the presence of columnar de- 
fects. The kinetic part of the hamiltonian with IMF is a convection-diffusion operator; 
operators of this type occur as well in models of driven population dynamics g , || , 
[Q. Due to the convection term, parity is violated and the wave function is expo- 
nentially magnified along a given direction; as a consequence, it can keep a localized 
character provided its decay, generated by the disordered potential, wins over the 
IMF amplification. 

In the zero-IMF reference hamiltonian the decay, if disorder does not violate parity, 
is characterized by the inverse of a single localization length on both sides. With 
convection, when the decay is exactly compensated by the IMF magnification factor, 
the wave function becomes delocalized. The spectrum at the critical point bifurcates, 
and is distributed along a d=l curve in the complex plane. This scenario was proved 
to be valid with various potential distributions (gaussian , box, Cauchy); it has been 
confirmed for lattice models with disordered hopping § , § , , ||] , |9| , jlO) . 
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In the hcrmitean case, one can turn the eigenvalue problem into a Langevin equa- 
tion, and determine the spectral properties pd| | from the steady state distribution of 
the Langevin problem. 

Here we generalize this approach to complex- valued stochastic equations, as needed 
when dealing with non-hermitean hamiltonians. Rather than starting with a gaussian 
^-correlated potential, which would lead to a Langevin equation, we consider a piece- 
wise constant potential V(x), a model originally studied by Benderskii and Pastur 
[ fl2| in the hermitean case. The reasons for our choice will become clear in the sequel. 
If, e.g., V(x) — Vo, V\ with Vq = and V\ < 0, V(x) describes equally engineered 
trapping defects in an homogeneous medium. 

We assume that free segments and trapping segments are independent random 
variables, with fixed average length Iq, l\ respectively. Although the stochastic equa- 
tion is no longer of Langevin type, the probability distribution still satisfies a dif- 
ferential equation. This equation allows for an adiabatic approach, which would be 
impossible with white noise. We can thus reduce the phase space to a one-dimensional 
attractor C, strongly simplifying the search for steady state. 

The "dynamic variable" of the stochastic equation is the logarithmic derivative 
z{x) of the wave function. From the mean value < z(E) >, (E = E r + i ■ Ei) one 
extracts the inverse of the localization length, as well as a wavenumber. Generically 
we find that the localization length decreases when E departs from the real axis, and 
increases with E r . 

This behavior implies that, in the presence of an IMF, at a critical value of E r the 
system delocalizes; the contour lines of < z r (E) > determine then the mobility edge 
in the complex energy plane. Three different regimes (low, high and intermediate E r ) 
will be analyzed in the next three sections. 

2 Langevin formulation of the quantum eigenvalue 
problem 

The eigenvalue problem Hip = E ■ ip for a d=l Hamiltonian Ho = + V(%) with 

disordered potential V(x) translates into a Langevin equation for the logarithmic 
derivative z of the wave function ip: 

^ = -{z* + E) + V{x); z = t- (1) 
ax ip 

As is well known, one can limit the analysis to real-valued wave functions. The 
deterministic part of the equation is invariant under space inversion (z(— x) — —z(x)) 
; hence, if the statistical properties of V(x) are also invariant, the backward and 
forward ^-evolution of z are equivalent, i.e. they approach the same steady state. 
One defines the Lyapunov exponent j{E): 

,. < lnr(x) > 

x 

r 2 (x) = ii 2 (x)+ip' 2 {x) 
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where the average is over V(x); 'y(E) is always positive for disordered potentials with 
rapidly decaying correlations 113]. The following identity is shown to hold: 



-/(E) =< z > 



(2) 



The Lyapunov exponent is an index for the exponential growth rate of the wave 
function (ip(x) w exp[±j ■ x] x — > +00): here, by Oseledec's theorem, the plus sign 
always occurs, but for a single case, corresponding to the physical (square-summable) 
solution |D| . Consistently, one has indeed < z > > and a non-even steady state 
distribution P(z). 

Let us now examine how this picture is changed upon adding an external (constant) 
imaginary magnetic field (IMF). In the hamiltonian, this is mimicked by: H = —(-^: — 
a) 2 + V(x), implying z — > z — a in Eq. [I], which now is to be studied for complex 
z and E ( z = z r + i ■ Zi, E = E r + i ■ Ei). The space-inversion invariance is 
broken by the convection term a ■ 4- and different exponential growth rates are found 
at x — > +00 and x — > —00. On qualitative grounds, the asymptotic behavior is 
ip(±L) Kt exp[(± ■ a— < z r (E) >) • L] (L >> 1) and derealization occurs when 
< z r (E) > -a = 0. 

If one gauges away the IMF from the Schrodinger equation, Eq. [l] is restored, but 
the boundary conditions are changed: e.g., in the case of a double-sided problem, from 
ijj(—L) — a_, ip(L) — a+, one has ip(—L) = exp(a ■ L) ■ a_, ip(L) = exp(—a ■ L) ■ a.+ . 

We are not going to study the spectrum, but rather the vanishing of the "effective" 
Lyapunov exponent < z r (E) > ~a. In order to do that it is sufficient to study the 
steady state distribution P(z) = P(z r , zi) of Eq. |l| with a — 0, obtain < z(E) > and 
from it the contour lines of < z r (E) >. Notice that from the analysis a new object 
emerges, < zi(E) >, a sort of wavenumber of delocalized solutions. 

In the search for steady state, we adopt an adiabatic approach. This is possible 
with our model potential, due to a property of Eq. |l| with V(x) = 0. Let us discuss 
this property first. 

Equation [l] has two time-independent solutions (respectively stable and unstable 
critical point). If we denote with e = e r + i ■ one of the square roots of E, the 
critical points are z s = e; — i ■ e r and z u = — + i • e r ,( the suffixes s, u meaning stable 
and unstable respectively, with e; > 0). It is easily verified that the Equation can be 
explicitly solved, for arbitrary initial conditions. After some algebra one obtains the 
announced property: the solutions satisfy an exact relation of the form: 



where j r and 7* and Xq are constants. 

As a model potential, we take V(x) two-valued, piecewise constant, V(x) = 
(Vq, Vi). The lengths of the Vi intervals are two independent random variables with 
distributions Qi(x) — jr ■ exp(—ci ■ x), (k — ^-). 

In order to write the analog of Equation |l|, we first introduce an independent 
stochastic variable s(x), which assumes the values 0, I with distributions Qo(x), Qi(x). 



\z(x)\ 2 



cosh(2B) - cos(2A) 
cosh(2B) + cos(2A) 
e r ■ (x - x ) + 7 r 
e, • (x - x ) + 7, 



(3) 



B 
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The probability Pq(x), Pi(x) of s(x) satisfies the rate equation: 



l£ = ^-A-Co-fl, (4) 

— — = c • Po - Cl • Pi, 

dx 



since s(x) stays on the 1-th channel an average "time"^-. 
The equation we looked for is then: 



dz 
dx 



(z 2 + E) + V + s(x)-(V 1 -V ), (5) 



and, in components: 



-z z r + zf -E r + V + s(x) • (Vi - Vq) (6) 



dz r 
dx 

d Z i r, 77. 



dx 

As is well known, from Eqjl], if V(x) is gaussian, one gets a Fokker-Plank equation 
for the probability. Equation ^ is no longer of Langevin type, but one can show 
that the probability distribution still satisfies a differential equation. The details of 
the derivation can be found in Ref. and in the original paper Jl2[ . In terms of 
the probability P(s,z;x) = Pq(z, x), Pi (z, x); (Pq(z;x) = P(s = 0,z;x), Pi(z\x) = 
P(s = 1, z; x) one has: 

dP ° - I d d p(°hp -Lr- P r p 

—r~ — (--, — i* r ' — 1—^1 )Po + ci • Jt\ — co • -ro 

ax dz r azi 

^ = (--^fW-A^) )Fi _ Ci .p 1 + Co . Po 

dx dz r azi 

= -zf + zf - E r + V (7) 

F^> = -zf + zf -Er + Vi 

ff = F^ = -2-z r - Zl - El 

We finally state our adiabatic approximation. We assume that the coefficients Co 
and Ci are very small with respect to the characteristic "frequencies" of the two 
deterministic evolutions z^°\x) and z^(x), respectively associated with E^ = E r — 
Vo+i ■ Ei and E^ = E r — V\ + i ■ Ei. Such frequencies, as one easily realizes, are 
the square roots e (0) , e (1) of E^ and E {1 \ 



3 Low energy case 

Let us examine the low E r case. In both deterministic systems we assume |e r | < |e$|. 
This regime is out of the physical region since it implies E r < Vi, (I = 0, 1); in spite 
of that, its analysis will be found to be useful in the sequel. 
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Upon averaging Eq. || over the variable B, which is fast varying with respect to 
A when |e r | < |ej| , the factor e°^(2B)+cos(2A) reduces to unity, and we obtain: 

1^(0(2 = | e (0|2, (/ = 0,1) (8) 

Each channel has a (metastable) attractor of circular shape, and the stochastic mo- 
tion reduces to hops between such concentric orbits. The single channel deterministic 
dynamics is illustrated in Figs.l and 2, where we show the vector field and the prob- 
ability distribution for Eq. || when s(x) is kept constant. 

We now proceed to determine the steady state of Eq. |7|. We first eliminate the 
Zi-dependence from Fr -* and Fr by means of Eq. ||, then average over Zi, and 
obtain: 

-^-(-2.z r 2 + 2.(ej 0) ) 2 )Po = CfP x -c -P (9) 

az r 

A(_ 2 . z 2 + 2- {ef^Pi = -a-Px + co-Po- 

Here Pq = Po(z r ) and Pi = Pi(z r ), are averaged over Zi . The sum of the two equations 
gives: 

(2 ■ zl - 2 ■ (ef ) 2 )Po + (2 • zl - 2 • (ef ) ) 2 )P 1 = K, (10) 

where K is a constant. 

It is readily verified that only with K = the positivity of the P's is preserved. 
When \z r \ is large enough the flow is in the negative z r direction on both channels; it 
turns positive for \z r \ < \ef'\ (I — 0, 1) respectively. 

Let us assume, without loss of generality, < < e . A simple inspection 
of the hopping between channels shows that the particle gets trapped in the interval 
ej < z r < e- ^; the probability is zero outside. Notice that upon including an 
indeterminacy in the radii of the two circular attractors we would get a non zero 
probability flow. 

Apart from the normalization factor H, the solutions have the form: 



POM = IT 1 L 

VI - (q 'r)\ 



P^Zr) = H- 1 m lP(z r ) 



_ (0) _ (1) 

*w - i\ z -^\r ■ [\ z -^\r (id 



Mi 



4 ■ e« 



The function p has an integrable power law singularity at the Z-th channel's stable 
point, and a zero at the opposite channel's stable point. If we denote with po and p± 
the integrals of Po(z r ) and Pt(z r ), they must satisfy the global equilibrium condition: 

co -po = ci -pi. (12) 
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We determined the solution in the region Co + c\ = 1 under this condition. The 
power exponent of the distribution is a function of the ratio j 2 -, i.e. it is shaped 
by the average "residence times" over the two channels. The mean value < z r > is 

increasing with the ratio p ~ -fjj, as shown in Fig. 5. Since larger ratios imply larger 

values of \Vg — V±\, this simply means that the Lyapunov exponent increases with 
disorder, as it should be. 

Any boundary condition z r {L) — ~z~Z with ~z~[ lying outside the interval e^p < z r < 
cannot be fulfilled: the density of states must then be zero in this regime, as 
anticipated. 



4 High energy case 

We analyze now the regime E r >> 0, which corresponds to the condition | | < |e r | 
for both channels. To simplify the notation, unless when strictly necessary, in the 
next formulas we drop the channel index (I). Contrary to the former case, in Eq. |^ 
the variable A is now fast varying with respect to B. The average over A gives: 

N 2 - n N 2 + H 2 (13) 

|e| 2 |e r ■ Zi — €i ■ z r \ ' 

from which: 

e| 2 = |e r • Zi - ti ■ z r \. (14) 

In a single channel, the metastable attractor is a couple of parallel lines: one through 
the stable, the other through the unstable critical point. Furthermore, the lines are 
orthogonal to the segment connecting such points. It appears that, in going from the 
low to the high energy limit, the attractor undergoes dilatation in the direction of 
such lines. This effect can already be seen in Fig. 2, referring to the low energy regime. 
An indeterminacy arises, about the line currently occupied by the "particle" : 

Zi = — -Zr±— (15) 

The single channel situation is illustrated in Figs. 3 and 4, where we display the vector 
field and the probability distribution for Eq. [| in the deterministic case. 

The stochastic motion involves the union of two couples of parallel lines, with 
distinct slopes. We proceed with our strategy, based on averaging Eq. fj] over Zj. 
Upon writing Ft (I — 0, 1) by means of Eq. [l^ we have: 

F r = -zl + (^ • z r f ±2-e l -z r + (e. ( ) 2 , (16) 

where the dependence on E and Vi is written in terms of e and the following approx- 
imation is made: ^ 

Zi = — ■ z r ± (e r + - — ) « — • z r ± e r . (17) 
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The force F^{z r ) {I = 0,1) acquires a fluctuating term, originated in the indetermi- 
nacy of Eq. [H| the 
force has the form: 



nacy of Eq. |15|; the term is of the order 2 • ef ■ z r . The dynamics associated with this 



d7 

= -z 2 r + (- ■ z r f + (af + 2-e t -z r - V (x), (18) 
ax e r 

where rj(x) is a white noise. 

In summary, the mentioned indeterminacy adds to the stochasticity of the hopping 
between channels I = 0, 1; the process is then described by the Fokker-Plank equation: 



dPo 
dx 


= A(_a(°) 
az r 


+ A jB (o)) Po + Cl . Pl _ Co .p 

az r 


dFi 

dx 


dz r ^ 


+ -^-B W )P 1 - ci • Pi + c ■ P 

az r 


ACQ 


= -{at-Zr) 2 


+ (4°) 2 


B (D 


= D-^r- 


[(^) 2 + (£ (i) ) 2 ] 




= i-(^) 2 





(19) 



where D is the white noise coefficient and £'Q is a regularizing parameter (the criteria 
used in fixing this parameter are illustrated in the caption of Fig. 6). 

We evaluated numerically, by a perturbative procedure, the steady state solution 
of Eq. II9L To simplify things, we assumed Co = cy = c < 1. When the channels 
are uncoupled (c = 0), the problem can be explicitly solved along the lines of the 
hermitean case (see, e.g., Ref. [11 ). The two equations reduce to the form: M^Pi = 



-Ji(z r ) — 0, where, with obvious notation, M^ is the second order operator acting 
on P t . 

The solutions correspond to a constant probability current J\ and here look like: 

Pi(zr) = [ r dz'exp[Ui(z r )-Ut(z')] (20) 

2 

1 ~1 



Bi{Zr) J- 



U t (z) = (l/D)[^p±-artg(- 



Upon substituting the series expansion Pi = J2k ° k ' m > = 0> 1) m ^1- @ 
one obtains an inhomogeneous equation for P^ k+1 \ where P, is the source term. 
The operator My' is easily inverted, with two integrals over z r . The integrals are 
performed numerically. In order to check for convergence, the area difference between 
successive approximants P" = J2k ° k ' * s computed. In the cases exibited here, 
we put c = 0.1, and find good convergence for n > 20. 

From the probability distribution Pi(z) we then get < z r >=< z r >o + < z r >x, 
with obvious notation. 

In Fig. 6 we plot the result < z r > as a function of E r . The asymptotic behavior in 
the hermitean case is known: < z r >« (E r )~ 2 , (Ei = 0) E r >> 1 juj|. Our estimate, 
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with Ei = 10. and E r in the range: 50. < E r < 120. , is of a much slower decay; 
further work is needed to clarify this point. 

When E r is fixed, < z r > increases with \Ei\, as shown in Fig. 7; we get < z r >« 
\Ei\P, f3 w 2 for large enough \Ei\. 

The value of the wavenumber < Zi(E) >, emerging in the delocalized regime, can 
be roughly estimated from < z r > and by choosing, in Eq. [15], the line through the 
stable point. 

It is easily verified that < Zi(E) > is an odd function of Ei while < z r (E) > is 
even; complex conjugate eigenvalues carry opposite wavenumbers. In the hermitean 
case the two eigenvalues merge , and the wavenumber, being a byproduct of broken 
parity, is zero for standing waves. 



5 Intermediate energies 

So far we studied two extreme cases, but intermediate energies in strong disorder can 
as well be treated; let us consider the following regime: Vj. << E r « Vo- 

The deterministic solutions satisfy Eqs. || and |l4] respectively, i.e. channel (0) is 
in a "low energy" regime, and channel (1) in a "high energy" one. The metastable 
attractor C is the union of a circle and of a couple of parallel lines. The system for 
Po and Pi is an hybrid between Eqs. || and [H]. 

The first order operator of Eq. ^ is coupled with the second order one of Eq. [l9]: 
this makes an iterative procedure hardly convergent. Based on results from direct 
numerical integration of Eq. ^|, we assign an indeterminacy S to the circular orbit: 

jfP = l + f. (21) 

An estimate for 5 can be extracted jD| from the deterministic flow (see the thick- 
ness of the ring in Fig. 2). The steady state equation then involves two second order 
operators: 

Ci • Pi - Co ■ P 

-ci • Pi + c • P (22) 



We proceed again by an iterative perturbation scheme. At uncoupled channels (c = 0) 
the "low energy" equation coincides now with the well-known one (jOl), holding in 
the hermitean case. Its solution carries a non zero probability current: in the present 
context, it comes from the spread of the circular attractor. The "high energy" integral 
was already exhibited in the previous section. 

In Fig. 9 we show the result, the surface < z r (E) >. The high energy regime surface 
is reported for comparison in Fig. 8; its contour lines in the E plane are also shown in 
Fig. 10. Qualitative agreement is found with previous numerical and analytical results. 
An estimate of < 2, (E) > can be derived from the adiabatic equations E5| and 0. 



-f [F (o)_-fnp o 

V° = (6-\e°\ 2 ) 2 



s 



6 Conclusions 



We studied nonhermitean delocalization in the presence of disordered potentials of 
Kronig-Penney type, where the barriers have random length. We translated the hamil- 
tonian eigenvalue problem into a stochastic equation for the logarithmic derivative of 
the wave function (z = ^-). From the steady state distribution of z we determined 
the Lyapunov exponent. 

Our approach, valid under suitable adiabatic conditions, describes the steady state 
in terms of the real part of z. The imaginary part of < z >, an effective wavenumber, 
can be recovered by means of the adiabatic formulae [sUnj 

We analyzed two regimes: a)intermediate energies (V\ « E r << Vq); b) high 
energies(Vi , Vq << E r ). The behavior of the system in other regimes is currently 
under investigation fl5|| . 

In case (a), the particle hops in the z plane between a strip and a circle, the strip 
carrying a probability current. The circle, which acts as a trap for the particle, is 
absent in case (b), where the 'particle' hops between two strips, with different slopes. 
Two currents coalesce in this case. 

We obtained that < z r > is always decreasing with E r : at higher energies the 
wave functions have larger localization length. 

At fixed < E r >, < z r > increases with \Ei\: this means that dissipation (associ- 
ated with Ei) enhances the localization. The effective wavenumber < zi > increases 
with E r , as one would expect with normal dispersion. 

The curve < z r (E) > —a = 0, a being the convection coefficient, determines the 
mobility edge in the complex energy plane. 

We finally add some comments on the density of states (DOS). Most analytical 
results on the DOS in the presence of an IMF, possibly with the single exception of the 
semiclassical analysis by Silvestrov , are concerned with the discrete case. Disorder 
averaging of the resolvent operator can indeed be performed in various discrete models, 
thus obtaining the DOS in explicit form @],§,|§,@. 

The Langevin approach gives the DOS for the hermitean hamiltonian in the con- 
tinuum. As a preliminary step in that procedure, one writes the definition of the DOS 
in terms of the z variable: this simply amounts to requiring that z(x, E) must fulfill 
the boundary condition: z(L,E) = z~Z- This condition is all is needed to determine 
the DOS, as long as the mapping from E to z(L, E) is invertible, which is precisely 
the case in the hermitean problem. 

Whether anything similar holds in the complex case is, as far as we know, an open 
question. 
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Fig. 1: Vector field (F r ,Fi) in the z plane for the Equation |B| at fixed s(x) (deter- 
ministic), in the low energy regime. The variables z r , Zi are in arbitrary units 
of (length)^ 1 . The superimposed circle corresponds to the metastable attractor 
discussed in the text (see Eq. ||). The two critical points characterized by the 
condition F r = 0, F% = belong to the circle. 
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Fig. 2: Numerical integration of the Equation |^ with s(x) = in the low energy 
regime: contour lines of the probability distribution at an intermediate time. The 
variables z r , zi are in arbitrary units of (length)^ 1 . The distribution is obtained 
from 2000 random initial conditions. The inner circle is superimposed: it cor- 
responds to the metastable attractor, as calculated from Eq. ||, with the values 
E r — Vo = — 56, Ei = —1.35 used in the numerical integration. The variables E r , Ei 
are in arbitrary units of (length)^ 2 . The unstable critical point is on the left of the 
ring; the peak at the stable point can be seen on the right. The distribution deviates 
from a purely circular behavior; it is elongated in the aequator direction, if one takes 
the critical points as the poles. In fact going from the low to the high energy regime 
(see Eq.|l5]) the attractor is converted into a couple of lines. These lines are through 
the poles, and parallel to the aequator. 
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Zr 

Fig.3:Vector field (F r , Ej) of Equation [| with fixed s(x) in the high energy regime. 
The variables z r , Zi are in arbitrary units of (length)^ 1 . The two superimposed par- 
allel lines through the critical points (F r = 0, Fi — 0) correspond to the metastable 
attractor (see Eq. |l5|). The lines are orthogonal to the segment connecting the 
critical points. 
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Fig. 4: Numerical integration of Eq. ^| with s(x) = 1 in the high energy regime 
(random initial conditions): contour lines of the probability distribution at an 
intermediate time. The variables z r , Zi are in arbitrary units of (length)^ 1 . The 
energy values are E r — Vi = 10.35, E.- L = 7.5, in arbitrary units of (length)^ 2 . The 
unstable critical point is in the upper left. The metastable attractor (two parallel 
lines) has been superimposed; with respect to Fig. 3, we are here at lower energies (the 
critical points are closer to the axis z r ). One can notice, in particular around the un- 
stable point, that the distribution is slightly elongated in the direction of the attractor. 
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Fig. 5: Low energy regime; from the solution (Eq. fyl) of the steady state Equation 
n e <°) 

|9| for the probability: log- log plot of < z r > versus p — -jjj. Higher p means higher 
disorder (see text): the localization length decreases with p. 
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Fig. 6: High energy regime: < z r > versus E r at fixed Ei, (Ei = 10). The variable 
z r is in arbitrary units of (length) -1 , and E in arbitrary units of (length) -2 . The 
origin of the energy axis corresponds to E r = | • (Vm) + V(i))', V(o) ~ V(i) = 5. From 
numerical integration of the probability Equation|l9| at steady state, with D = 0.1. 
The regularization parameter is: = e\. With this choice £^ has the same order 
of magnitude of z r at the stationary points: z r = ±e| . 
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Fig. 7: High energy regime: log-log plot of < z r > versus E{ at fixed 
E r : E r — V(i) = 60., E r — V(o) = 30. Same Equation and parameters as in 
Fig.6. 
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120 

Fig. 8: High energy regime:global behavior of < z r > versus E; V(p) — V(i) = 30. The 
variable z r is in arbitrary units of (length)^ 1 , and E in arbitrary units of (length)^ 2 . 
From numerical integration of Equation [l9] for the steady state probability. Same 
parameters as in Figs. 6,7. 
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Fig. 9: Intermediate energy regime, from numerical integration of the Equation |22| for 
the steady state probability: < z r > versus E,D = 0.1, 5 = 25., Vrm — Vm = 30. The 
variable z r is in arbitrary units of (length) -1 , and E in arbitrary units of (length) -2 . 
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Fig.lO:High energy regime. Contour lines of the surface < z r (E) > plotted in Fig. 8. 
E ri Ei are in arbitrary units of (length)^ 2 . Each line identifies a mobility edge, for 
a given convection parameter a; a increases from bottom right to top left, in a range 
of values from 0.1 to 2.95, with step 0.15. 



20 



